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O , We found the parameterization of the unintegrated gluon distribution from the best de- 

scription of the LHC data on the inclusive spectra of hadrons produced in pp collisions at 
the mid-rapidity region and small transverse momenta. It is different from the one obtained 
within perturbative QCD only at low intrinsic transverse momenta kt- The application 
of this distribution to analysis of the e — p DIS allows us to get the results which do not 
contradict the HI and ZEUS data on the structure functions at low x. So, the connection 
between the soft processes at LHC and low-x physics at HERA is found. 
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As is well known, hard processes involving incoming protons, such as deep-inelastic lepton- 
proton scattering (DIS), are described using the scale-dependent parton density functions. Usu- 

p^ ally, these quantities are calculated as a function of the Bjorken variable x and the square of the 

four-momentum transfer q^ — —Q^ within the framework of popular coUinear QCD factoriza- 
tion based on the DGLAP evolution equations [I]. However, for semi- inclusive processes (such 
as inclusive jet production in DIS, electroweak boson production [2], etc.) at high energies it is 

^\ more appropriate to use the parton distributions unintegrated over the transverse momentum 

^ kt in the framework of fcj-factorization QCD approach |3 F] The /c(-factorization formalism is 

based on the BFKL |6j or CCFM [7j evolution equations and provides solid theoretical grounds 
for the effects of initial gluon radiation and intrinsic parton transverse momentum kt- The 
theoretical analysis of the unintegrated quark q{x,kt) distribution (u.q.d.) and gluon g{x,kt) 
distribution (u.g.d.) can be found, for example, in [SJ- [24]. According to 11 , the u.g.d. 



g{x,kt) at fixed Q^ has the very interesting behavior at small x < 0.01, it increases very fast 
starting from almost zero values of kt and decreases smoothly at large kt- In contrast, the 
u.q.d. q{x, kt) is almost constant in the whole region of kt up to kt '^ 100 GeV/c and much 
smaller than g{x,kt). These distributions were obtained using the so-called KMR formalism 
within the leading order (LO) and next-to-leading order of QCD (NLO) at large Q^ from the 
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^See, for example, reviews [4 5 for more information. 



MPI@LHC 2011 



known (DGLAP-evolved [T]) parton densities determined from the global data analysis. The 
unintegrated parton distributions were succefully applied to analyze the DIS data at low x and 



the number of processes studied at LHC (see, for example, 13 - 23 . However, at small values 
of Q^ the nonperturbative effects should be included to evaluate these distributions. The non- 
perturbative effects can arise from the complex structure of the QCD vacuum. For example, 
within the instanton approach fl2] the very fast increase of the unintegrated gluon distribution 
function at < fcf < 0.5 GeV/c and Q^ = 1 (GeV/c)^ is also shown. These results stimulated us 
to assume, that the unintegrated gluon distribution in the proton can be included by analyzing 
also the soft hadron production in pp collisions. In this paper we analyze the inclusive spectra of 
the hadrons produced in pp collisions at LHC energies in the mid-rapidity region including the 
possible creation of soft gluons in the proton. We estimate the u.g.d. at low intrinsic transverse 
momenta kt < 1.5 — 1.6 GeV/c and its parameters extracted from the best description of the 
LHC data at low transverse momenta pt of the produced hadrons. We also show that our u.g.d. 
similar to the u.g.d. obtained in l9)fTo] at large kt and different from it at low kt- 



2 Inclusive spectra of hadrons in pp collisions 

2.1 Unintegrated gluon distributions 

As was mentioned above, the unintegrated gluon densities g{x, k^) can be generally described 
by the BFKL [6] evolution equation, where the leading ln(l/a;) contributions are summed. 
The conventional integrated gluon distribution g(x,Q^) can be approximately related to the 



unintegrated one by 13 



g{x, Q') - g{x, Ql) + / dk^tgi^x, k^t) (1) 



where Qq is a starting nonzero value of Q^. An appropriate description valid for small and 
large x is given by the CCFM [T] evolution equation that results in the u.g.d. g{x, kf,q^) as a 
function of a;,fc^ and the additional scale q (see details in [4J[lO] and references therein). The 
another u.g.d. g{x, k^) satisfies the DGLAP-type evolution equation with respect to k^ [2 
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A simple parameterization of the u.g.d. was obtained, for example, within the color-dipole 
approach in [8l[9] on the assumption of a saturation of the gluon density at low Q^ which 
succcefuUy described both inclusive and diffractive e — p scattering. The u.g.d. xg(x, fcf , Qq) is 
given by [9 10 



xg{x, kt,Qo) = -—^ — TTTT^t^ '^^P (^^0(2^)^?) 5 ^0 = 7^ ( — ) , (2) 

where a^ = 29.12 mb, a^ = 0.2, Qq = 1 GeV/c, A = 0.277 and xq = 4.1- 10"^ The form 
for xg{x, kt,Qo) given by Eq.([2| was obtained in [9] within the model of the {qq) dipole 24 25 
on the assumption of the saturation effect for the gluon density, e.g., the dipole- nucleon cross 
section a^^.^ is assumed to be a constant at low Q^ (and low x too). It corresponds to the 
Gaussian form for a^^pflh) as a function of the impact parameter b. In fact, this form could 
be more complicated. In this paper we study this point and try to find the parameterization 
for xg{x, kt, Qo), which is related to a-y^]\j(h), from the best description of the inclusive spectra 
of charge hadrons produced in pp collisions at LHC energies and mid-rapidity region. 
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2.2 Quark-gluon string model (QGSM) including gluons 

As is well known, the soft hadron production in pp collisions at not large transfer can be analyzed 
within the soft QCD models, namely, the quark-gluon string model (QGSM) [26]- 28 or the 
dual parton model (DPM) [29]. The cut n-pomeron graphs calculated within these models 
result in a reasonable contribution at small but nonzero rapidities. However, it has been shown 
recently 



30 



that there are some difficulties in using the QGSM to analyze the inclusive spectra 
in pp collisions in the mid-rapidity region and at the initial energies above the ISR one. It is due 
to the Abramovskiy-Gribov-Kancheli cutting rules (AGK) [3l] at mid-rapidity (y — 0), when 
only one-pomeron Mueller-Kancheli diagrams contribute to the inclusive spectrum ph{y — 0,pt). 
To overcome these difficulties it was assumed in [30] that there are soft gluons in the proton, 
which are split into qq pairs and should vanish at the zero intrinsic transverse momentum 
{kt ~ 0). The total spectrum Phiy — 0,Pt) was split into two parts, the quark contribution 
Pq{y ~ 0,p() and the gluon one and their energy dependence was calculated [30] 



p{Pt) = Pq{x = 0,pt) + Pg{x = 0,pt) -gis/sn) <j)q{0,Pt) + {g{s/sQ - and) (pg{0,Pt) ■ 



The following parameterization for 0^(0, pt) and (j)g{0,pt) was found 

4>q{0,Pt) = Aqexp{-bqPt) 

4>g{0,Pt) = Ag^te^i-bgPt), 



30 



(3) 



(4) 



where sq = IGeV^, g = 21m&, A = 0.12. The parameters are fixed from the fit to the data on the 
Pt distribution of charged particles at y = [30]: Ag = 4.78 ±0.16 (GeV/c)-^, 6, = 7.24 ±0.11 



(Gev/c)-i and Ag == 1.42 ± 0.05 (GeV/c)-^; bg = 3.46 ± 0.02 (GeV/c)-i. Figure 1 illustrates 
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Figure 1: The inclusive spectrum of the charged hadron as a function of pt (GeV/c) in the 
central rapidity region (y = 0) at y/s =7 TeV at pt < 1.6 GeV/c compared with the CMS 32 
which are very close to the ATLAS data [33] . 

the best fit of the inclusive spectrum of charged hadrons produced in pp collisions at ^/s —7 
TeV and the central rapidity region at the hadron transverse momenta pt < 1.6 Gev/c; the 
solid line corresponds to the quark contribution pq, the dashed line is the gluon contribution 
Pg, and the dotted curve is the sum of these contributions ph given by Eq.ra. 
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2.3 Modified unintegrated gluon distributions 

We calculated the gluon contribution 4>g{Q,Pt) entering into Eq.W as the cut graph of the 
one-pomeron exchange in the gluon-gluon interaction (Fig. 2, right) using the splitting of the 
gluons into the qq pair. Then the calculation was made in a way similar to the calculation of 
the sea quark contribution to the inclusive spectrum within the QGSM, see Eqs.(4,5) in 30 at 
n = 2. 
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Figure 2: The one-pomeron exchange graph between two gluons in the elastic pp scattering 
(left) and the corresponding cut graph (right). 



Pg{x±,Pht) ^ F{x+,Pht)F{X-,Pht) , 



(5) 



where 



F{X±,Pt;Pht) ^ dxi(fkltfq,q{xi,klt)Gg(^q)^hi—,Pht-kt)] , (6) 

Here Gqiq\^h{z, kt) — zDqtq\^fi{z, fcf ), Dqtq\^i^{z, kt) is the fragmentation function (FF) of the 
quark (antiquark) to a hadron h, z — x±/xi,kt — Pht ~ kt, x± — Q.b{^Jx'^ + Xj ± x),xt = 
2yJ(rn\ + p1)/ s. The distribution of sea quarks (antiquark) fqq is related to the splitting 
function Vg^qq of gluons to qq by 



fq..q{z,kt) 



g{zi,kt,Qo)Vg^qq{ 



Zl 



(7) 



where g{zi, ku, Qo) is the u.g.d.. The gluon splitting function Vg^qq was calculated within the 
Born approximation. 

Calculating the diagram of Fig. 2 (right) by the use of Eqs.(l5|7]) for the gluon contribution pg 



we took the FF to charged hadrons, pions, kaons, and pp pairs obtained within the QGSM 34 
From the best description of pg{x ~ 0,Pht, see its parameterization given by Eq.Q, we found 
the form for the xg(x, kt, Qo) which was fitted in the following form: 



xg{x, kt, Qo) = , /''L , ^1(1 - x)"^ X 
AiT^asiQo) 

{Rl{x)kl + C2{Ro{x)kt)'') exp {-Ro{x)kt - d{Ro{x)ktf) , 



(8) 
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The coefScient Ci was found from the following normalization: 



Jo 



(9) 



and the parameters 

a = 0.7; C2 ^ 2.3; A = 0.22; bg = 12; d = 0.2; C3 = 0.3295 

were found from the best fit of the LHC data on the inclusive spectrum of charged hadrons 
produced in pp collisions and in the mid-rapidity region at pt <1.6 GeV/c, see the dashed line 
in Fig. 1 and Eq.Q. 
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Figure 3: The unintegratcd giuon distribution xg{x,kt,QQ)/Co as a function of kt ai x — xq 
and Qo = l.GeV/c. The dashed curve corresponds to the original GBW [9{[To], Eq.([2]), and the 
solid line is the modified u.g.d. given by Eq.®. 



Figure 3 presents the modified u.g.d. obtained by calculating the cut one-pomeron graph of 

[9 10 as a function of the transverse gluon momentum kt. 

(the solid line in Fig. 3) is 



Fig. 2 and the original GBW u.g.d 

Here Cq — 3cro/(47r^as(Qo))- One can see that the modified u.g.d 



different from the original GBW u.g.d. [9 10 at kt < 1.5 GeV/c and coincides with it at larger 
kt- This is due to the sizable contribution of pg (Eqs.(3J4|) to the inclusive spectrum p{pt) of 
charged hadrons produced in pp collisions at LHC energies and in the mid-rapidity region, see 
the dashed line in Fig.l. 
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3 Proton longitudinal structure function 

Within the fct-factorization approach, the proton longitudinal structure function F]^(x,Q^) 



calculated in the leading order of QCD can be presented in the following form 13 



Fl(x, 



dz 

z 



dk 



flavour{f) 



(10) 



r,m}, 



k1) are the so 



where e^ is the charge of the quarks of flavor /, the functions C|^(a;/2; 

called hard structure functions of the off-shell gluons with virtuality fc^ which correspond to the 
quark-box diagram for the photon-giuon interactions |13j . In the present note we calculated 
Fl{x, Q^) at the fixed value of the missing mass W — 276 GeV using the parameterization for 
u.g.d. at fixed Qo given by Eq.Q and Eq.dsl. The results of our calculations are presented in 
Fig. 4 in comparison with the HI data 35 36 
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Figure 4: The longitudinal structure function i^L(/x|j) at W =276 GeV and /i^ = K -Q^, where 
X = 127 l37|. The HI data are taken from f35)|36]. 



We find that the modified u.g.d. allows us also to describe the proton Fl Note that in order 
to take into account the NLO corrections (which are important at low Q^) in our numerical 
calculations we apply the method proposed in Il4j. Following 14 37 , we use the shifted value 
of the rcnormalization scale /i|j, — K Q^, 



37 



this shifted 



where K = 127. As is was shown in 
scale in the DGLAP approach at LO approximation leads to the results which are very close 
to the NLO ones. 
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4 Conclusion 

We fitted the experimental data on the inclusive spectra of charged particles produced in the 
central pp collisions at energies larger than the ISR starting with the sum of the quark con- 
tribution pq and the gluon contribution pg (see Eqs.(3J4|). The parameters of this fit do not 



depend on the initial energy in that energy interval. Assuming creation of soft gluons in the 
proton at low transverse momenta kt and calculating the cut one-pomeron graph between two 
gluons in colliding protons we found the form for the unintegrated gluon distribution (modified 
u.g.d) as a function of x and kt at the fixed value of Qg. The parameters of this u.g.d. were 
found from the best description of the LHC data on the inclusive spectra of the charged hadrons 
produced in the mid-rapidity pp collisions at low pt. It was shown that the modified u.g.d. is 
different from the original GBW u.g.d. obtained in [9,10 at kt < 1.6 GeV/c and it coincides 



with the GBW u.g.d. at kt > 1.6GeV/c. It was also shown that the modified u.g.d. allows us 
to describe the longitudinal structure function Fl{Q^) at the fixed missing mass W better than 
the original GBW u.g.d. Therefore, some link between soft processes at the LHC and low- a; 
physics at HERA is found. 
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